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Harmonic maps on domains with piecewise 
Lipschitz continuous metrics 

Haigang Li* Changyou Wang ^ 



Abstract 

For a bounded domain Q. equipped with a piecewise Lipschitz continuous 
Riemannian metric g, we consider harmonic map from (Q, g) to a compact Rie- 
mannian manifold (N, h) without boundary. We generalize the notion of 

stationary harmonic map and prove the partial regularity We also discuss the 
global Lipschitz and piecewise C l a -regularity of harmonic maps from (fl, g) to 
manifolds that support convex distance functions. 

1 Introduction 



Throughout this paper, we assume that Q. is a bounded domain in R", separated by a 
C u -hypersurface Y into two subdomains Q + and Q.~, namely, Q. = Q. + U £l~ U T, and g 
is a piecewise Lipschitz metric on Q. that is g e C 01 (£2 + ) n C 0,1 (Qr) but discontinuous 
at any x e Y. For example, Q. = B x c R" is the unit ball, r = B x n {x = (x', 0) 6 R"}, 
• ' and 

(go xeB+ = {x n >0}nB 1 , 

g(x) = i 

[kgo xeB- = {x n <0}nB l , 

where go = dx 1 is the standard metric on R" and k 1) is a positive constant. 

Let (N, h) M, k be a /-dimensional, smooth compact Riemannian manifold with- 
out boundary, isometrically embedded in the Euclidean space R*. 

Motivated by the recent studies on elliptic systems in domains consisting of com- 
posite materials (see Li-Nirenberg IfTTlO and the homogenization theory in calculus of 
variations (see Avellaneda-Lin HI and Lin-Yan lfT8l0 . we are interested in the regular- 
ity issue of stationary harmonic maps from (Q, g) to (N, h). 

In order to describe the problem, let's first recall some notations. Throughout this 
paper, we use the Einstein convention for summation. For the metric g = g t j dx 1 dx j , 
let (g' ; ) denote the inverse matrix of (gij), and dv g = yfgdx = Jdet (gy) dx denote the 
volume form of g. For 1 < p < +oo, define the Sobolev space 

W lp (Q, N) = jw : Q. -» R* u(x) e a.e. xeQ, E p (u, g) = J \Vu\ p g dv g < +oo| , 
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where 



._ l2 u . du du v 
5 <7JQ cuc ; - 



is the L 2 -energy density of w with respect to g, and <•, •) denotes the inner product in 
R*. Denote W l2 (Q,N) by H^QN). 

Now let's recall the concept of stationary harmonic maps. 

Definition 1.1. A map u e H 1 ^, N) is called a (weakly) harmonic map, if it is a 
critical point of E 2 (-, g), i.e., u satisfies 

A g u + A(u)(Vu, Vu) g = (1.1) 

in the sense of distributions. Here 

is the Laplace-Beltrami operator on (f2, g), A(-)(-, ■) is the second fundamental form of 
(N, h) ^ R\ and 

A(u)(Vk, V«) g = g y A(«) ^— , — , 

Definition 1.2. A (weakly) harmonic map u e H l (Q.,N) is called a stationary har- 
monic map, if, in additions, it is a critical point of E 2 (-,g) with respect to suitable 
domain variations: 



d 
dt 



|V«f dvg = 0, with u\x) = u(F t (x)), (1.2) 



where F(*,x) := F t (x) e C 1 ([-6, 6], C 1 

is a C 1 family of differmorphisms for 

some small 6 > satisfying 

F (jc) = jc Vjc 6 Q, 

F f (x) = x V(x, f) 6 x [-5, £], (1.3) 
F,(?F)c?F Vf £[-£,<$]. 

It is readily seen that any minimizing harmonic map from (Q, g) to (iV, /z) is a 
stationary harmonic map. It is also easy to see from Definition 1.2 that a stationary 
harmonic map on (f2, g) is a stationary harmonic map on (Or,g) and hence satisfies 
an energy monotonicity inequality on Or, since g e C°' l (Qr). We will show in §2 that 
a stationary harmonic map on (Q, g) also satisfies an energy monotonicity inequality 



in Q under the condition (|1.4I) below. 

The first result is concerned with both the (partial) Lipschitz regularity and (par- 
tial) piecewise C ^-regularity of stationary harmonic maps. In this context, we are 
able to extend the well-known partial regularity theorem of stationary harmonic maps 
on domains with smooth metrics, due to Helein 021, Evans flU, Bethuel Q. More 
precisely, we have 
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Theorem 1.1. Let u e H l (Q.,N) be a stationary harmonic map on (£l,g). If, in addi- 
tions, g satisfies the following jump condition on Y for n > 3 |jj : for any x e F, there 
exists a positive constant k(x) ± 1 such that 




lim g(y) = k(x) lim g(y), 



(1.4) 



then there exists a closed set IcU with H n 2 (2) = 0, such that for some < a < 1, 



(i) u e Lip loc (ft \ E, AO, (ii) u e Cj£((Q + U T) \ 2, AO n c£(tflT U T) \ S, AO. 



We would like to remark that when the dimension n = 2, since the energy mono- 
tonicity inequality automatically holds for i^-maps, Theorem l 1 . 1 I holds for any weakly 
harmonic map from domains of piecewise C^-metrics, i.e., any weakly harmonic map 
on domains with the above piecewise Lipschitz continuous metrics is both Lipschitz 
continuous and piecewise C l,a for some < a < 1. 

Through the example constructed by Riviere [|T9l , we know that weakly harmonic 
maps on domains with smooth metrics may not enjoy partial regularity properties in 
dimensions n > 3. Here we consider weakly harmonic maps on domains with piece- 
wise Lipschitz continuous metrics into any Riemannian manifold (N,h), on which 
djj(-,p) is convex. Such Riemannian manifolds Af include those with non-positive 
sectional curvatures, and geodesic convex balls in any Riemannian manifold. In par- 
ticular, we extend the classical regularity theorems on harmonic maps on domains 
with smooth metrics, due to Eells-Sampson [8] and Hildebrandt-Kaul-Widman |fT3l , 
and prove 

Theorem 1.2. Let g be the same as in Theorem \l.l\ Assume that on the universal cover 
(N,h) of {N, the square of distance function d-(-, p) is convex for any p e N. If 
u 6 H l (Q., N) is a weakly harmonic map, then for some < a < 1, 



The idea to prove Theorem 1.1 is motivated by Evans (H and Bethuel 0. How- 
ever, there are several new difficulties that we have to overcome. The first difficulty 
is to establish an almost energy monotonicity inequality for stationary harmonic maps 
in Q, which is achieved by observing that an exact monotonicity inequality holds at 
any x e T, see §2 below. The second one is to establish a Hodge decomposition in 
L P (B, W), for any 1 < p < +oo, on a ball B(= 5, (0)) equipped with certain piecewise 
continuous metrics g, in order to adapt the argument by Bethuel [2J. More precisely, 
we will show that the following elliptic equation on B: 



This condition is needed for both energy monotonicity inequalities for u in dimensions n > 3 and 
the piecewise C l,a -regularity of u. 

2 Here the covering map n : — > is a Riemannian submersion from (A^, h) to (A^, h). 



(i)«eLip loc (Q,A0, 



(ii) u e c!' a (Q + ur,N)n cjf (*r u r, AO. 




= on dB 
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enjoys the W^-estimate: for any 1 < p < +00, 



Vv < C 

[ \\LP(B) 



f 



LP(B) 



provided that (a*/) 6 cf#*) n C^B 6 ^) for some 8 > is uniformly elliptic, and is 

discontinuous on dB + \ B s , where B 6 = 6 B : dist(jc, 55) < 

This fact follows from a recent Theorem by Byun-Wang y), see §3 below. The 
third one is to employ the moving frame method to establish a decay estimate in suit- 
able Morrey spaces under a smallness condition, which is similar to [fT4l . To obtain 
Lipschitz and piecewise C x - a -regularity, we compare the harmonic map system with 
an elliptic system with piecewise constant coefficients and extend the hole-filling ar- 
gument by Giaquinta-Hildebrandt IfTOl . 

The paper is organized as follows. In §2, we derive an almost energy monotonicity 
inequality. In §3, we show the global W 1,p (1 < p < 00) estimate for elliptic systems 
with certain piecewise continuous coefficients, and a Hodge decomposition theorem. 
In §4, we adapt the moving frame method, due to Helein [fT2l and Bethuel B2]|, to 
establish an £-Holder continuity. In §5, we establish both Lipschitz and piecewise 
Ci,a regularity for Holder continuous harmonic maps. In §6, we consider harmonic 
maps into manifolds supporting convex distance functions and prove Theorem 1 1 .21 
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NSF grant 1000115. 



2 Energy monotonicity inequality 

This section is devoted to the derivation of energy monotonicity inequalities for sta- 
tionary harmonic maps from (Q., g) to (N, h). More precisely, we have 

Theorem 2.1. Under the same assumption as in Theorem \l.l\ there exist C > and 
r > depending only on Y and g such that ifu £ W l ' 2 (Q, N) is a stationary harmonic 
map on (Q, g), then for any xo e Q., there holds 

\Vu\ 2 dv g <e Cr r 2 - n f \Vu\ 2 dv„ (2.1) 

JB r (xo) 

for allO < s < r < min{r , dist(x , dQ.)}. 

Since the metric g e C 01 ^*), it is well-known that there are K > and r > 
such that (12.11) holds for any xq 6 and < s < r < min{r , dist(xo, dQ.*)}, see lfT2l . 
In particular, (12.11) holds for any xq 6 Q. \ Y r ° and < s < r < min{r ,dist(xo,5Q)}, 
where P = |x e : dist(x, T) < r } is the r -neighborhood of Y. We will see that to 
show (12.11 ) for xo e F" , it suffices to consider the Case Xq € 1 . 
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It follows from the assumption on Y and g, there exists r > such that for 
any xo 6 T there exists a C^^differmorphism O : B{ — » B ri (xo), where ri = 
min{r , dist(x , 5f2)}, such that 

1 (Ti) = T n B n (jcb), where H = {* 6 ^ : x„ = 0}. 

Define w"(x) = i/(O (x)) and g(x) = (<b ) t (g)(x) for jceBi. Then it is readily seen that 

(i) gis piecewise C 01 , with the discontinuous set H, and satisfies (|1.4I) on T\ [E 

(ii) If u : (B n (x ), g) — » (Af, /z) is a stationary harmonic map, so does u : (Bu 'g) — » 
(N, h). 

Thus we may assume that Q. = By, g is a piecewise C 01 -metric which satisfies 



(|1.4I) on the set of discontinuity r l5 and u : (B u g) — » (Af, /z) is a stationary harmonic 
map. It suffices to establish (|2.1I) in Bi. We first derive a stationarity identity for u. 

Proposition 2.2. Le? u e W 1,2 (fii,A0 oe a stationary harmonic map on (B u g). Then 

X (<<■ £>* - *H - X, v^k <£. £>* m 

holds for allY = (Y 1 ,-- - , Y n ~\ Y n ) e Cj(fli, R") satisfying 



Y n (x) 



> /or > 

= for x" = (2.3) 
< /or jc" < 0, 



dY^ " ^ 
where Yf = and div Y = > . 

(=i 

Proof. Let F satisfy (12.31) . it is easy to see that there exists 6 > such that F t (x) = 
x + tY(x), t 6 [-5,5], is a family of differmorphisms from B\ to B\ satisfying the 
condition (1.3). Hence 



°4 

dt 



^ JT |V(u(F,(*))g Jv g = ^[_ () ( |V(«(F,(x))i;dVj + |V(«(F,(x))g</v,). 



d® k d¥, 

In fact, since (<&o)*(g)ij(x) = gki(^o(x))— — (x)— — (x), ( 11.41 ) implies that for any x e H 

OX, OXj 

lim (G> o )»g()0 = *(*o(*)) lim (O )^O). 

yen + ,_y— >.v yeu. ,y — >jc 



5 



d 
dt 
d 
dt 



For t e [-6, S], set G t = F t l . Direct calculations yield 

r iv( M (F,w)i^ Vg 

yfg(^g i} (x)(^, ^r>(* + tY(.x))(S u + tYf)(6u + tY l j) dx 
'=oJs± oy k dyi 1 

= [ rS8 ij {^-,^-){S k1 Y l j + 6 li Yl)dx 
J B ± dx k dx, ] 

+ f t\ {g i \G t {x))^^{GM)JG t {x)){^-,^-)dx 
J B ± dt\t=o v ' dxi oxj 

C n,du du^, •• du du s \ 
Jg+5x/t^^^ ^ ^dx;' dx) 

1 ^ 



where we have used 



JG t (x) = -divF, 
G t (x) = -Y(x), 



t=0 {g ij (G t (x)) V^W)) = -^{yfg^Y 
This completes the proof. 



□ 



Proposition 2.3. Let u e W li2 (Bi,N) be a stationary harmonic map on (Bi,g). Then 
there exists C > such that 



(i) for any x u = (x Q , 6 fii \ fi, there exists < R < min{|, \x^\), such that 



"I 

JB r (jfi) 



\Vu\ 2 g dv g < e CR R 2 ~ n \ \Vu\ 2 g dv g , < r < R < R . (2.4) 



(ii) for any x° £ B\ n Y\, there holds 

r 2 ~" f IVwgrfVg < e CR R 2 ~ n f \Vu\ 2 e dv 
JbjjP) " Jb, 



g ^ g , 0<r<R<- (2.5) 



In particular, for any x° e Bi, there holds 



i 



\Vu\]dv g < e CK R 



CRr,2 



i 



\Vu\ 2 g dv g , 0<r<R<-. (2.6) 



Proof (i) By choosing Y e C~(5|,R M ) or Y e C^(B^,R n ), we have that u is a sta- 
tionary harmonic map on ( B\, g) and (B^,g). Thus the monotonicity inequality (12.41) 
is standard, see |fT2l. 
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(ii) For simplicity, consider x° = (0',0). For e > and < r < \, let Y e {x) = 
xt] £ (x), where rj £ (x) = rj £ (\x\) e C™(Bi) satisfies 



< T] £ < 1; r] e (s) = 1 for < s < r - e; rj e (s) = for s > r; rj' € < 0; \tj' £ \ < -. 



Then 



(n)/ = ^(W) + ^(W) 



XX 



Substituting Y £ into the right hand side of (|2.2I) . and using 



dx k 



we have 



Substituting (12.71) into the left hand side of (12.21) . we obtain 

= (2-n) f IV^WV?^- f |V M |Jjc|^(jc)VgdJc 



+ 



C ~ a,du d u s xk x J , , n r- , 

2 S ,7 <t-,t->tt^VI^- 



Set the piecewise constant metric g by 

g(x',x n ) 



lim g(y) if jc" > 

y-»0, y>o 



lim e(v) if x n < 0. 

ly->0, >>"<0 



Then we have 



\g(x) - g(x)\ < C\x\, VxeB,. 
It follows from (11.41) that we can assume 



g(x) 



\go if ^ > 
\kg ifx n <0, 



for some positive constant k £ 1. Thus we can estimate 
,du du v , 

/B, 



(2.7) 



(2.8) 



(2.9) 



(2.10) 



f „ , du du 4 x^x- 7 . _ , 
2^<— ,— >— -rf € {x)^gdx 
J Bl dxi dx k \x\ 

„ C _ij,du du x k x J . _ , _ f , ; ; -ih,du du x k x j v 

J Bl <9x, dx^ |x| J B] dx,- dx* |x| 



(2.11) 
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Since 



_y du dtl Jt*X> = (\x\lff if X" > 

8 { dxi'dx k } \x\ ~ {X) - \{\x\\ d f r \ 2 ifx n <0, 



and h(x) > for x e B\, we have 

h= f h(xW £ (\x\)^dx<0. (2.12) 
Jb, 



For 7/ e , by (12101) we have 



• f |Vw| 2 dv g <Cr \ 

*J By *J B r 



<CY I \Vu\ z dv„<Cr I |V a | 2 Jv„. (2.13) 



First substituting (12.121 ) and (12.131 ) into (12. lib , and then plugging the resulting (12.1 II ) 
into (12.91 ), and finally combining (12.91 ) and (12.81 ) with (12.21) , we obtain, after sending e 
to zero, 

(2-n) f \Vu\ 2 g dv g + r f \Vu\ 2 g Vs dH n ~ l > -Cr f \Vu\ 2 g dv g . 

JB r JdB r JB r 

This implies 



|VH|^ Vg |>0, 



which clearly yields (12.51) . 

To show (12.61) . it suffices to consider the case 

x° 6 B m \ T u \B R (x°) n > and |5 r (jc°) n 5^1 > 0. 

For simplicity, assume x° e 5 r We divide it into two cases: 
(i) d(x°,r 1 ) = \x° n \>lR: 



4 

2- 



• If R > r > ^R, then it is easy to see 

~ n \ \Vu\ 2 g dv g < A n - 2 R 2 ' n \ \!vtf g dv g . 

If < r < \R(< d(x°, TO), we have Br(x°) c B~ so that (Q implies 
r 2 "" f |V M grfv g < e CR (j) f |V«|^Vg < e CR R 2 - n [ \Vu\ 2 g dv g . 



(ii) d( J c°,r 1 ) = |^|<i/?: 



• If 7? > r > then 



r 2 ~" f \Vu\ 2 g dv g < 4 n ~ 2 R 2 -" f \Vu\ 2 g dv g . 

JB r (x°) " JB R (x°) 
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• If < r < d(x°, TO = \x°\ < \R, then by setting x = (x° v • • • , x° n _ v 0) we have 
B r (x°) c B^(x°) c B M (x°) c Br c 
so that ((231) yields 

r i i » hi. ,n . ,, _i i.»„i i i » ii|,,n * 



f \Vu\ 2 g dv g < \x Q f- n f m 2 g dv 

<\r n r 

(T r 



< 2"- 2 (2|4 , |) 2 "" I ^ |V M |^ 

7 x ■* 



<e CR R 2 ~ n \Vu\ 2 g dv g . 

If d(jc°,ri)(= < r < \R, then we have 

B r (x°) c fl 2 r(*°) c c B R (x°), 

so that (1231) yields 



- f \Vu\ 2 g dv g <2"- 2 (2r) 2 -" f \Vu\ 2 g dv 

{-)"" f 



,2-n 

in— 2 „CR | 



< 2"--/ A (-| I ^ |V«|^v g 



Jb r (x0) 



<e CR R 2 ~ n | \Vu\]dv g . 



>B R (x») 

Therefore (12.61) is proven. □ 

3 W 1 ^-estimate for elliptic equations with certain piece- 
wise continuous coefficients 

In this section, we will show the global Vl^-estimate for elliptic equations with cer- 
tain piecewise continuous coefficients, for 1 < p < +oo. As a corollary, we will estab- 
lish the Hodge decomposition Theorem 3.2 for certain piecewise continuous metrics 
g, which is a key ingredient to prove Theorem [Tj] and may also have its own interest. 

For a ball B = B r (0) c R", denote B e = {x e B : dist(x, dB) < e} for e > 0. 
Let (ciij(x))i<jj< n be bounded measurable, uniformly elliptic on B, i.e., there exists 
< A < A < +oo such that 

m 2 < «//(-v)^> ^ m\ a.e. xeB,^e R". (3.1) 
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Theorem 3.1. Assume satisfies d3.il) . and there exists e > such that (a, 7 ) G 
C (2F) n C (5 e ) and is discontinuous on dB + \ B e . For 1 < p < +oo, /g? / e L P (B, R n ). 
Then there exists a unique weak solution v G W^' P (B, R") to 

r Z|-(^|-) = Zf inB, 

' i 

on dB, 



u = 



(3.2) 



and 



l|Vv| 



L"(B) 



< c 



L"(B) 



(3.3) 



/or some C > depending only on p and (a,;). 

Proof. By our assumption, it is easy to verify that for any 5 > 0, there exists i? = 
R{5) > such that the coefficient function (a, ; ) satisfies the (5, i?)-vanishing of codi- 
mension 1 conditions (2.5) and (2.6) of Byun-Wang [3] page 2652. In fact, we have a 
stronger property: 



lim max 

40 X0 =( x ' Q ^)eB 



Qij(x , X ) aij(X(\, xf 1 ) 



= 0. 



L~(B r ((y ,^))) 

Thus Theorem 13.11 follows by direct applications of [3] Theorem 2.2, page 2653. □ 

As an immediate consequence of Theorem 13 . 1 L we have the following Hodge de- 
composition on B equipped with certain piecewise continuous metrics g. 

Theorem 3.2. Let g be a piecewise continuous metric on B such that g e C {b*} H 

c(B s )for some 5 > 0, and is discontinuous on dB + \ B s . Then for any 1 < p < +oo, 

F = (F u • • • , F„) G LP(B, W), there exist G G W^' P (B) and H G LP(B, IT) swc/j f/raf 



1 

F = VG + H, = div^tf (:= — — ( Vfs' 7 #,)) in 5, 



(3.4) 



and ?/zere exists C = C(p, n, g) > smc/z ?/za? 

||VG|| L;)(B) + ||//|| LP(B) <C||F|| i;)(B) . (3.5) 

Proof. Set a iy = -^gg' j on 5 for 1 < i,j < n. It is easy to verify that (a, ; ) satisfies 
the conditions of Theorem 13.11 Thus Theorem 13.11 yields that there exists a unique 
solution G G W l p {B) to 

{A(VIf i7 f) = A(VI^ F y) ^ 5, 

I G = on 55, 



(3.6) 



and 



l|VG|| 



L»(B) 



< C 



V(B) 



< c 



L"(B) 



Set H = F - VG. Then we have 



div s H 



1 8 



dG 



dxj 



0on£, 



and 



\H\\v (Bi ) < \\F\\ mBl ) + \WG\\ LP(B) <C\\F\\ LP(B) 



This completes the proof. 



□ 
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4 Holder continuity 

In this section, we will prove that any stationary harmonic map on (Bi,g), with a 
piecewise Lipschitz continuous metric g e C°' l (B* U Ti), is Holder continuous under 

a smallness condition of J |Vm|^ dv g . The idea is based on suitable modifications 

of the original argument by Bethuel [2] (see also Ishizuka-Wang lfT4l ). thanks to the 
energy monotonicity inequality and the Hodge decomposition theorem established in 
previous sections. More precisely, we have 

Theorem 4.1. There exist eo > and ao e (0, 1) depending only on n,g such that if 
the metric g e C°' l (Bf U Y[) satisfies the condition H1.4H on Y[, and u e W (B\, N) is 
a stationary harmonic map on (B\,g) satisfying 

-" f \Vu\ 2 g dv g <e 2 (4.1) 

JB rQ (x ) 

for some xq e B\_ and < r Q < j, then u e C a °(B[o(xo), N) and 

< C(r , e ). (4.2) 



M 



(x )) 

Proof of Theorem \4J\ The proof is based on suitable modifications of [2] and lfT4l . 
First, observe that if x = (xL, x'£) 6 5*, it follows from the monotonicity inequality 
(|2.6I) that we may assume (14.11) holds for some < r < \x'l\. Then the e -regularity 
theorem by Bethuel [21 (see |[l4l for domains with C ' 1 metrics) implies that for some 
< ao < 1, u e C a{) {Bn(x )) and (14.21 ) holds. Hence it suffices to consider the case 

2 

x = (x' Q , 0) e r i . By translation and scaling, we may assume x$ = (0, 0) and proceed 
as follows. 

Step 1. As in |[2l ffT2l IfPfl . assume that there exists an orthonormal frame on u*TN . 
For < 6 < \ to be determined later, let {e a } l a=l c W 1 ' 2 (B 2 e, R k ) be a Coulomb gauge 
orthonormal frame of u*TN : 

Bio 

| div g ((Ve a , e p )) = in B 2e (1 < a,B < I), 
\ij B2 Ve a \]dv g < C f B JVu? g dv g . (4J) 

For 1 < a < I, consider (V ((w - u 2 g)r]) , e a ), where u 2 e = T u is the average of u on 
B 2e , and rj e Cq(Bi) satisfies 

2 

< 77 < 1; 77=lin5 e ; rj = outside Bi e ; |V//| < -. 

4 

Let go be the standard metric on R". We define a new metric ~g on B 2e by letting 
J(x) = r](x)g(x) + (1 - r](x))g (x), x 6 ^2^. 
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Then it is easy to see that 

gE^onfi 9 jEg outside Bi_ g , and g e C(B± e ) n C(B 2e \ Bi_ ). 

In particular, ~g satisfies the condition of Theorem 13.21 Hence, by Theorem 13 .21 we 
have that for 1 < p < -, there exist <b a e w!. ,p (B 2e ) and \b a e L p (B 2 e) such that 



(V ((u - u 2 s)n) . e<*) = V^ ff + ifr a , divg(^) = in B 2e , 

< l|Vw|| 



(4.4) 



Since u satisfies the harmonic map equation (11.11) . we have 

6i\ g «Vw, e a )) = g' J ViU{Vje a , e p ) e p in B y . (4.5) 

Thus we obtain 

A g <p a = g' J ViU{Vje a , e p ) e p in B e . (4.6) 



Set (p a = (p)/ + (p)/ , where (p a l) solves 



and <^ 2) solves 



[aX 1) = mB e , 
\(p ( a=(p a ondBo, 



^g<Pa ] = g ij Viu(Vje a , e p ) e B in B e , 
[4 2) = ondB e . 



(4.7) 



(4.8) 



Step 2. Estimation of <f$\ It is well-known (cf. 0H) that 0^ e C ao (fi e ) for some 
o-o e (0, 1), and for any < r < | 

KLm,^ ^ f IV^Y^^C^"" f |V M |^x, (4.9) 

and 



*J B T ti 



where M p ' p (-) denotes the Morrey space: 



, VO < t < 1, (4.10) 

Mi>-P(B,) 



M p > p (E) :={f:E^R: 11/11'' = sup {r p ~ n f \f\ p dx\ < +oo), E c R\ 

B r (.v)cR" ^ JB r (x)nE ) 

Step 3. Estimation of (p^: First, denote by "H^R") the Hardy space on R" and 
BMO(£) the BMO space on E for any open set E c R". By (4.13) of flHl page 
435, for p = — ^-j- > n, there exists h e W^ P \B ), with ||V/j|| L ,/ (Be) = 1, such that 

< C f (V(pf,Vh) g dv g . 

J Bf) 
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Hence by the equation (14.81) . (14.41) . and the duality between and BMO, we have 

M 2) L ( s 9 )- C f <S8 ij ^iU^je a ,e p )){e p h)dx 

J Bg 



£ II V^^< V 7^' e /?)lli 2 (B«)ll V (^)lli 2 (B«) MbM0(B 2 «) 

< l|VM|| L 2 (B2fl) ||V a || MP ,, (Bl) • 0?"*, 



= -c 

< C 



S29) 



(4.11) 



where we have used: 

(i) Since div ? «Ve Q ,, e p )) = in 5 and h e W ' p (B 9 ), we have y/gg' J (Vje a , ep))V i{eph) € 
<H\R n ) and 



V&"< V a , ^>v ( -M) wi(Rn) < c || Vim*- ^>|| i2(Be) ||v ; -M)|| L2(B( 



l— — 

(ii) Since /?' > n, the Sobolev embedding implies h e C ? and 



so that 



W(efih)\\i? m < liv^ii^^pn^^) + \\vh\\ D , m e-r-^ < ce~r\ 

(iii) By Poincare inequality, it holds 

MBMO(B 2fl ) - c H Vm IImp/'(b,)- 

Putting the estimates of <pa and 0^ 2) together, we obtain 

1 

[ijey-' 1 ^ \V<f> a \ p dxJ < C [r a ° + r'-T'eo] Wu\\ MP , nBl) , VO < r < 1. (4.12) 
Step 4. Estimation of i// a : Since diVg^Q,) = on B 2e , we have 

I \^a\^dVj= J ((lf/ a + V(f> a ),lf/ a )-dVg 
<J Boa *J Biq 



= I ((V((u-u 2 e)T]),e a ),if/ a )gdVg 
Jb 2 b 

= - I (u - u 2e )r]{Ve a , ip a )-gdvg 
Jb 1b 



< 



gg tJ ^ i e a i// J a [(u - u 2e )T]] BMO 



e a\\L-(B 2l j) [(" - M 2e)77] BM o 

£ l|V«|| L 2 (B2fl) ||lA Q ,|| Z ,2 (B2e) ||VM|| MW(Sl) , 
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where we have used the fact 

[(u - u 2 o)t]] bmo < C Mbmo(b 28 ) < C ||Vw|| 
This, combined with Holder's inequality, implies 



< Ce \\Vu\\ 



(4.13) 



Step 5. Decay estimation of Va: Putting (14.121) and (14.131) together, we have that for 
some < a < 1, 



\(r9y- n j \Vu\ p) j" < C(e + r a ° + tH^^Vm 



Mi'P(Bi) 



(4.14) 



holds for any < t < 1 and < 6 < |. Now we claim that for some a £ (0, 1), it 
holds 



\\Vu\\ MP , HBr) < C(e () + r a ° + r 1 -^,,) IIV 



Mll^fp.n-p^j) , V0 < T < 1. 



(4.15) 



To show (14.151) . let B s (y) c Br. We divide it into three cases: 

(a)yeBrC\B ± and s < \y n \. As remarked in the begin of proof, we have that for some 
< or < 1, 



Jb, 



(v) 



iv»l"| <c[- 

s 



< C — 

Vi 



< c 



\y n r n f |V« 

(2\f\r n f 

xz/ JSrfy',0) 



|V M |' 



(since \y"\ < -) 



< C(6 + r a ° + T l -J>e Wu\\ M , P(Bl) (by 
(b) y £ Bi n B ± and s > \y"\. Then we have B s (y) c B w+s (y',0) c B 2s (y',0). Hence 



f-n f 
Jb, 



(v) 



Vu\ p \ <2—\(2s) p 



JB 2s (y',0) 



< C(e + T ao +t 



J ~-'eo)\\Vu\\ MP , P{Bl) (by (El). 



(c) y e Br_n T u i.e. y" = 0. Then it follows directly from (|4.14l) that 

i 

f |V«n P < c(e + r ao + T l -f>e ) ||V«|| MP , P(Bl) . 

Combining (a), (b) and (c) together and taking supremum over all B s (y) c Br, we 
obtain (14.151) . 
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It is now clear that by first choosing sufficiently small r and then sufficiently small 
6q, we have 

I|Vm]|mp.P(Bt) < — ]|Vw]| M ;),;.( Bl ) . 

Iterating this inequality finitely many time yields that there exists a\ e (0, 1) such that 
for any x e Bi and < r < \, it holds 



[ \Vu\ p dx<Cr^\\Vu\\ p MP ,, (B{) . 

JBAx) 



Ax) 

This immediately implies u e C ai {B\_). The proof is now completed. □ 

5 Lipschitz and piece wise C 1 ^-estimate 

In this section, we will first establish both Lipschitz and piecewise C 1 '" -regularity for 
stationary harmonic maps on domains with piecewise C 01 -metrics, under a smallness 
condition of energy. Then we will sketch a proof of Theorem ll.il 

Theorem 5.1. There exist e > and B e (0, 1) depending only on n,g such that if 
the metric g e C°' l (Bf U Ti) satisfies the condition < \1.4\) on Y[, and u e W (B\, N) is 
a stationary harmonic map on (B\,g) satisfying 



[ \Vu\ 2 g dv g <e 2 (5.1) 

JB r Jxo) 



for some x Q e B^ and < r Q < |, then u e C l ^°^Bn(x Q ) n i^Af), and u e 
C°' 1 (b^(x ),N). 

Proof. The proof is based on both the hole filling argument and freezing coefficient 
method. It is divided into two steps. 

Step 1. u e C a (Biro{xty), N) for any < a < 1. To see this, recall Theorem 14. 1 1 implies 

4 

that there exists < «o < I sucn that u e C a °(B Tro(xo)) and for any y e Bt^(xo), it 

J 8 8 

holds 



s 2 ~ 



" f |Vw| 2 dx < C I-) 2 " r 2 -" f \Vu\ 2 dx, < s < r < ~, (5.2) 
J BAy) X - r ' J BAy) 8 



3Ay) x,/ -J BAy) 

and 



8 

For y 6 5 v, (x ) and < r < let v : B r (y) -> R* solve 

(A g v = in £ r (y) 
|v = w on dB r (y). 



osc Br(y) u < Cr ao , < r < ^. (5.3) 



(5.4) 
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Then by the maximum principle and (15.31) . we have 

osc Br0 ,)V < osc dBr(y) u < Cr ao . 

Moreover, since g e C 0,1 (B* U H), it is well-known (cf. IfTTII Theorem 1.1) that 
v e C 0A (B^(y), R k ) and v e C l *(Br(y) n IF, R*) for any < B < 1. 

Now multiplying both the equations (11.11) and (15.41) by (u - v) and subtracting each 
other and then integrating over B r (y), we obtain 

f |V(» - v)| 2 Jjc < f \Vu\ 2 \u -v\< r"- 2+3ao . 

JB r (y) JB r (y) 



Since 



we obtain 



I 



\Vv\ 2 dx<C 



Vv 



|V M | 2 < C( ||Vv|| 2 „ (Br W) r 2 + r 3 «°) < CV 3 T 



3gQ 

This, combined with Morrey's decay lemma, yields u e C 2 (Bir (xo)). Repeating this 



8 

-<ar 



argument, we can show that u e C (B3r (x )) for any < a < 1, and 

4 

r 2 ' n f \Vu\ 2 dx < Cr 2a , Vv e B*Axo), < r < -. (5.5) 

JB r (y) " 4 4 

Step 2. There exists < j8 < 1 such that u e C ia (b^(xo) n IF,/v). The proof is 
divided into two cases. 

Case I. x = (x^xf^) e B*. We may assume < r < \x^\ so that B ro (xo) c B ± . For 
B r (x) c S ro (jc ), let v : 5 r (jc) -> R* solve 

(A g v = in fl r (*), 
)v = w on dB r {x). 

Then by Step 1, we have that for any | < a < 1, 

f |V(m-v)| 2 Jjc< C f IVi/plM-vlJjc^Cr 3 ^"- 2 . (5.7) 
Moreover, since g e C^CB^Qto)), we have that for any < B < 1, v 6 C l ^{B^{x)) and 
f |Vv-(Vv) Bs( ,)| 2 d*<C(-) 2/ * f |V« - (V M ) BrW | 2 0<s< r -. (5.8) 
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Henceforth, we denote T f = — \ f dx. Combining (15.71) and (15.8lH we obtain 

Je \E\Je 

that for any < 6 < 1, 

\ \Vu - (Vu) Berix) \ 2 dx<l\\ |Vw - Vv| 2 dx + \ |Vv - (Vv) B9rW | 2 dx] 

< C\e 2fi f \Vu - (Vu) Br(x) \ 2 dx + 6' n r 3a - 2 ]. 



>B r (x) 

For ^=1 < B Q < 8, let < 6 Q < 1 be such that Cflf = ef °. Then we have 

f |Vu - (Vuh 8oAx) \ 2 dx < 6f° f \Vu - (Vu) BAx) \ 2 dx + Cr ia ~ 2 . (5.9) 

JB Sor (x) JbAx) 

Iterating (15.91 ) m-times, m > 1, yields 



f |v M -(v M ) BvW | 2 jx < (^) 2/3o f |v« - (V W ) Br( ,)| 2 <fc 

JB g m r (x) JB r (x) 

m 

+C(ff£r) 3a - 2 ^ ^»~( 3a - 2 » (5.10) 
< (^)M f |Vm - (V M ) BrW | 2 J* + Cr 3a - 2 ]. 

JB r (x) 

3a-2 

This clearly implies that Vu e C - (B ro (xo)). 

Case II. x = (x' ,0) e Ty. For simplicity, we assume x' () = 0. Define the piecewise 
constant metric g on B x by letting 



g(x) = 

Then we have 



[lim fT0 ^(0',0 xeB~. 



\g(x)-g(x)\<C\x\, x G B\, (5.11) 

Moreover, by suitable dilations and rotations of the coordinate system, (11.41) implies 
that there exists a positive constant k + 1 such that 

g(x) = (1 + (k - 1 )XB-(x))go, x 6 5i, 



whereas- is the characteristic function of B x . 

£0 

2 



For < r < ^, let v : B r (0) R k solve 



fA,v = mB r (0), (5i2) 
lv = w on cJ5 r (0). 



4 note that ( 15. 8t trivially holds for j < s < r. 
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Then we have 

osc fir(0) v < osc Br(0 )W < Cr a , \ \Vv\ 2 dx < C f \Vu\ 2 < Cr n ~ 2+2a . 

JB r (0) JB r (0) 

Multiplying (11.11) and (15.121) by (u - v) and integrating over B r (0), we obtain 
f \V(u-v)\ 2 dx 

Js r (0) 

< I g'\u-v)i(u-v)j^[gdx 

Jfi r (0) 

< C f \Vu\ 2 \u-v\dx+ f \^g ij - ^g ij \\vi\\(u-v)j\dx 

JB r (0) Jb,-(0) 

< Cosc B ,.(o)V f \Vu\ 2 dx + Cr 2 f |Vv| 2 + \ f \V(u-v)\ 2 dx 

< Cr"" 2 ^ + Cr n+a + \ f \V(u-v)\ 2 dx. 
This implies 

f |V(m-v)| 2 Jjc< Cr"- 2+3a . (5.13) 

Js,-(0) 

It is well-known that v 6 C°°(i?±(0)) for any < s < r. In fact, (15.121) is equivalent to: 

—((1 + (*i - VXb-)-^) = 0, in B r (0), (5.14) 

we conclude 

(i) Jj- satisfies the jump property on Ti : 

lim — (x',x n ) = lim -^(x',x n ), V(x\0) e^n fl,.(0). 

-v„io+ dx n jf„T0- obc„ 

(ii) V a v 6 C°(5 r (0)) for any multi-index or = (ac\, ■ ■ ■ , or„_i, 0). 

(iii) Vv 6 L°°(fl,(0)) for any < s < r, and 

HVv||^ (Br(0)) <Cr 2 "" f |V M | 2 . (5.15) 

2 JB r (Q) 

For / : 5,-(0) R*, set 

OX\ OX n ^\ ox n 

and let (Df) = \ Df dx denote the average of Df over B s (0). Then we have that 

V /s Jbao) 
for any < B < 1, 

-f |5v - (Dv) s \ 2 dx<C (-) -f |5m - (Du) r \ 2 dx, V0 < s < r. (5.17) 
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Combining (15.131) with (15.171) yields that for any < 6 < 1, 



f \Du - (Du) dr \ 2 dx < C0 2/? f \Du - (Du) r \ 2 dx + Cr'V 3 ^ 2 . (5 J 8) 
JB e M JBM 



As in Case I, iteration of (|5.18l) yields that for any < s < r, it holds 

■f \Du - (Du) s \ 2 dx < cl-V 2 f \Du- (Du) r \ 2 dx + Cs^ 2 . (5.19) 
This, combined with Case I, can imply that for any B r (x) c B ro (xo) and < s < r, 
f \Du - (Du) x ^ dx < c(-) f \Du - (Du)^\ 2 dx + Cs 3a ~ 2 , (5. 
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where (Du) xs denotes the average of Du over B s (x). It is readily seen that (15.201) yields 
Vw G C hJ ^ (B r_o(x ) n Bf) and u e C -\B[o(x )). This completes the proof. □ 

Now we sketch the proof of Theorem ll.il 

Proof of Theorem [7TT1 . Define the singular set 

2 = Left: limr 2 -" f \Vu\ 2 dx>e% 

{ r->0 JB r (x) 

Then by a covering argument we have H"~ 2 (Y) = (see, for example, Evans-Gariepy 
0). For any xq e ft \ S, there exists < r < dist(x , 5ft) such that 



|Vw| 2 JjC < 6q. 



rl~ n f 
Jb T0 ( X ) 

Hence by Theorem l2.ll Theorem l4.1[ and Theorem 15. 1[ we have 

u e C l ' a (Bm(x Q )nOF,N) and u e ^(b^xoXN), 
for some < a < 1 . In particular, we have 

lirnr 2 "" \ \Vu\ 2 dx = 0, Vx e B^(x ), 
so that B'_q (x ) n 2 = and hence 2 is closed. This completes the proof. □ 



6 Harmonic maps to manifolds supporting convex dis- 
tance functions 

In this section, we consider weakly harmonic maps u from (ft, g), with g the piecewise 
Lipschitz continuous metric as in Theorem 1.1, to (N, h), whose universal cover (N, h) 
supports a convex distance function square d~(-,p) for any p e N. We will establish 
both the global Lipschitz continuity and piecewise C 1 '" -regularity for such harmonic 
maps u. This can be viewed as a generalization of the well-known regularity theorem 
by Eells-Sampson [8 J and Hildebrandt-Kaul-Widman |fT3l . 

The crucial step is the following theorem on Holder continuity. 
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Theorem 6.1. Assume that the metric g is bounded measurable on Q., i.e. there exist 
two constants < A < A < +00 such that AI n < g(x) < AI„ for a.e. x e £1 Assume also 
that the universal cover (N, h) of (N, h) supports a convex distance function square 
d 2 -(-,p)forany p e N. If u e H l (Q,N) is a weakly harmonic map, then there exists 
a G (0, 1) such that u e C a (Q, N). 

Proof. Here we sketch a proof that is based on modifications of that by Lin lfT6l . 
Similar ideas have been used by Evans in his celebrated work |[6| and Caffarelli 01 
for quasilinear systems under smallness conditions. First, by lifting u : Q. — » N to 
a harmonic map 77 : Q, — > N, we may simply assume (N,h) = (N, h) and d 2 N {-,p) is 
convex on N for any p e N. 
We first claim that 

A g d 2 (u,p) > 0. (6.1) 
In fact, by the chain rule of harmonic maps (cf. Jost lfT5l ). we have 

A g d 2 (u,p) = V u d 2 (u,p)(A g u) + Vld 2 (u,p)(Vu,Vu) g . 

Since A g u ± T U N, V u d 2 (u,p) 6 T U N, the first term in the right hand side vanishes. By 
the convexity of d 2 N , the second term in the right hand side satisfies 

V 2 u d 2 (u,p)(Vu,Vu) g > 0. 

Since u 6 H l (£l, N), by suitably choosing p 6 N and applying Poincare inequality and 
Harnack's inequality, (16.11) implies u e LJ° c (£l,N). 

For a set £ c N, let diam A ?(£') denote the diameter of E with respect to the distance 
function d N (-, •)• For any ball B r (x) c Q, we want to show that u e C a (B>_(x)) for some 

< a < 1 . To do it, denote 

C r '■= diam/v (u(B r (x))) < +00. 

We may assume C r > (otherwise, u is constant on B r (x) and we are done). Now we 
want to show that there exists < 5q = 6 (N) < \ such that 

diam w (u(B Sor (x))) < ]-C r . (6.2) 

Since u r (y) = u(x + ry) : B\(0) — » N is a harmonic map (5i(0),,g r ), with g r (y) = 
g(x + ry), we may, for simplicity, assume x = and r = 2. For any < e < \, since 
w(5i) c N is a bounded set, there exists m = m(e) > 1 such that u{B\) is covered by m 
balls B 1 , • • • , B m of radius tCi. Now we have 

Claim: There exists a sufficiently small e > such that u(Bi) can be covered by at 
most (m — 1) balls among B , ■ ■ ■ , B'" . 

To see this, let x, 6 B { such that B' c B 2e ci(Pi), Pi = for 1 < i < m. Let 

1 < m' < m be the maximum number of points in {p, } j such that the distance between 
any two of them is at least ^C\. Thus B± Ci (pi), 1 < i < m', covers u(Bi). Then there 
exists i'o £ {1, ■ •• , m'} such that 

ic?<sup4(«(jc),/>fc)<C?, (6.3) 

4 jceB 2 
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and 



H n [u- 1 (B N (p i0 ,^C l ))nB 1 \> c () , 



(6.4) 



for some universal constant c > 0, where B N (p ia ,R) is the ball in N with center p i( 
and radius R. 
In fact, since 

Brclju-'^ipi,— d)j, 

m' I , \ 

2 ^" 1 (5 A '( A , —CO) n 5! j > 



we have 



Hence there exists io e {1, • • • , m'\ such that 



tf" (u\B N {p i0 , -Lei))) > co := —H n (B\). 
\ 16 / m' 



This implies (16.41) . By the triangle inequality, (16.31) also holds. 
Define 

/(x) := sup^(«(z),/7 io ) - d 2 N (u(x),p io ), xeBi. 

zefli 

It is clear that / > in B\, and (16.11) implies 

A g /<0, infli. 
By Moser's Harnack inequality, we have 

inf/>C f f >C \ f >C \ f 



finK-Wten.^Ci)) 
1 1 " 



> c 



mpd N (u,p io ) 



sup d 2 N (u,p io ) 



nf 1 ^,,- CO) 



1 
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^C(^-^C?]c :=« 



(6.5) 



for some universal constant O > 0. This implies 

sup d N (u(z),p k ) - sup d N (u(z),p io ) > 6 Q C\ = (1 - o )Ci. 



(6.6) 



Now we argue that the claim follows from (16.61) . For, otherwise, we would have 
that w(fii) n B 2eCl (Pj) + © for all 1 < j <m. Let z e #i be such that 

eCi + d N (u(z ), p io ) > sup d N (u(z),p k ). 
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Since u{B x ) c U'" = B 2£Cl (Pi), there exists p h e {p u ■■■ ,p m } such that u(z ) e B 2eCl (p h ). 
Since u(B i)nB 2e c l (Ph) ^ 0, there exists zi 6 B\_ such that u{zi) e B 2e Ci(Ph)- Therefore 
we have d N (u(z{), u(zo)) < 2eCi. Therefore we have 

sup d N (u(z),p k ) - sup d N (u(z),p k ) < 6d+ d N (u(z ), p io ) - d N (u(zi), p io ) 

< ed+ d N (u(zo), u{z\)) < 3eCu 



l 



this contradicts (16.61) if e > is chosen to be sufficiently small. 
From this claim, we have either 

(i) diamtf(K(fli )) < )-Ci. Then (Q holds with S = \, or 

(ii) diam N (u(B^)) > -Ci. Then we consider v(x) = : B { — » iV and have 

• v is a harmonic map on (2?i, gi), with the metric gi (x) = g(^rx). 

• ici < diam w (v(5i)) <Ci. 

• v(.Bi) is covered by at most (m - 1) balls B\, ■ ■ ■ , B m ~ x of radius eC\. 

Thus the claim is applicable to v so that u(Bi) = v{B\_) can be covered by at most 
(m - 2) balls among B 1 , ■ ■ ■ , B m '\ 

If diam^(v( J Bp) < -Ci, we are done. Otherwise, we can repeat the above argu- 
ment. It is clear that the process can at most be repeated m-times, and the process will 

not be stopped at step k < m unless diam N u(B 2 -k ) < -C\. Thus (16.21) is proven. 
It is readily seen that iteration of (16.21 ) implies Holder continuity. □ 

Proof of Theorem [772] First, by Theorem 6.1, and the argument from §4, we can show 
that for some < a < 1 , 



\Vu\ 2 dx < Cr n ~ 2+2a , MB r (x) c Q. 



lB r (x) 

Then we can follow the same proof of Theorem 5.1 to show that u e C 01 (fi) and 
u e C l ' a (Qr U T,N). □ 
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